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SUMMARY: The growth of pulmonary function between 5.5 and 25 yr 
of age was determined using 1,511 observations over time on 353 
subjects from a representative population sample of white non- 
Mexican-Aroericans in Tucson. There was an average of 8.8 yr of 
follow-up, with a maximum of 12. The method used was shown to be 
robust for span of follow-up from 3 to 12 yr (3 to 7 observations), 
and the results were verified by standard statistical methods. 
The standard error of the estimate decreased linearly with follow¬ 
up, indicating th eneed for longitudinal evaluation. Respiratory 
symptoms and diagnoses had the biggest negative impact on' growth' 
of lung function, using FVC, FEV1, Vmax50, and size-compensated : 
flows (Vmax50/FVC). Smoking had the next biggest negative impact. 
Smoking cessation was shown to have a positive impact on growth of 
pulmonary function. Using a second linear model to adjust for 
individual variability and the random variability over surveys, 
individual growth showed similar trends. Further negative impacts 
were due to parental smoking, especially as it interacts with active 
smoking and respiratory disease. Flows at end of follow-up (Vmax50, 
Vmax50/FVC) were more sensitive than FEV1 to the effects of 
concurrent disease and smoking, and more persistent effects of 
these factors in early adulthood. 
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Introduction 

I he study of the growth of pulmonary 
function in childhood and into early 
adulthood has stimulated many investi¬ 
gations. Although there have been vari¬ 
ous evaluations of this “growth” cross- 
seasonally, fewer studies have evaluated 
such growth through longitudinal mea¬ 
surements of pulmonary function (1-9). 
Many of these studies have been with 
only 2 observations over lime. For those 
studies with 3 or more points in time, the 
basic questions relate to how to repre¬ 
sent the generalized growth model, how 
growth is statistically dependent on ini¬ 
tial measures of function, how it is related 
to respiratory illnesses in childhood, and 
the effects of active smoking, passive 
smoking, and other environmental fac¬ 
tors. Basic questions are whether there 
are differences in growth between the 2 
sexes (10), and whether airways of differ¬ 
ent sizes grow at the same rate or are af¬ 
fected differently by internal and exter¬ 
nal factors (11). 

This report evaluates the longitudinal 
growth of pulmonary function in sub- 
)ccts who have been enrolled in the Tuc¬ 
son Epidemiological Study of Airways 
Obstructive Diseases since 1972 (12 )l It 
is a study of the development of spiro- 
metric function and the characteristics 
and robustness of a model of growth. The 
effects of sex, age, size, and of other re¬ 
spiratory faaors on growth are examined. 
The report attempts to evaluate the ef¬ 
fects of initial function, symptoms, dis¬ 
ease during childhood, the onset of 
smoking, parental symptoms and paren¬ 
tal smoking. 

Methods 

Subjects studied were those enrolled in the 
longitudinal epidemiologic study of white 
n on-Mexican-American households in TUc- 
SOn - The study has been described in detail 
elsewhere (12), 


Spirometric Testing Over Time 
Subjects 5.5 yT of age or older were tested 
spirometrically with a pneumotachygraph, 
using ATS/Snowbird criteria, as described pre¬ 
viously (13, 14). The measures of pulmonary 
function derived were FVC, FEV,, flow at 
50*7* of FVC (Ymax*,), and flow at 75^o 
of expired FVC (Vmax n ). Their percent of 
predicted FEV, values were in the range 
of 49 to 151^* at entry, using out previously 
published prediction formulks (14). 

Subjects were tested in Surveys 1 to 8 
(February 1972 through May 1984), with the 
exception of Survey 4. Thus, it was possible 
to have a maximum of 7 observations. We 
required at least 3 yr of follow-up for jub- 
jecu to be included in the data file. There were 
• total of 1,511 observations on 353 subjects 
5.5 to 15 yT of age at the time of their first 
satisfactory testing, and the oldest was 25 yr 
of age at the end of follow-up The maximum 
length of follow-up was 12 yr, and the aver¬ 
age was 8.8 yr; table 1 shows the MU obser¬ 
vations divided by the number of observations 
per person. The height distribution by age ap¬ 
peared normal for a pediatric population] 

The age range was divided into four 5-yr 
intervals to evaluate cohort effects; cohorts 
were divided into those in an age group who 
had their first spirometric test either during 
the first 3 yr of the study (1972 to 1975) or 
whose first test was performed during subse¬ 
quent surveys (1977 or later). Using multiple 


regression techniques, each pulmonary func¬ 
tion variable was regressed against age (years, 
decimalized)^ height (in inches), and the 2 to¬ 
gether, separately for those within each age 
group and by cohort. Differences were evalu¬ 
ated by multifactorial analyses of variance 
(ANOVA) as well, in which gender was a 
covariate. In these analyses, gender was not 
significant for Vmix w , and was inconsistent 
for FEV,. There were no significant differ 
cnees between cohorts in the specific regres¬ 
sions in the 4 age groups (5.5 to 10 yr v 11 to 
15, 16 to 20, and 21 to 25), for any of the 4 
measures (FVC, FEV,, Vmax*,, Vmax M ), as 
confirmed by ANOVA. The Vmax,, results 
were highly variable, and were not used fur*, 
then Thus, we studied all subjects with longi¬ 
tudinal data, regardless of when they entered 
the longitudinal study. 
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TABLE 1 

THE NUMBER Of OBSERVATHENS PER INDIVIDUAL SUBJECT 


2 3 4 5 6 7 Tout 


Ot>*#rv«tK>ns 

n 

56 

210 

440 

425 

240 

140 

1.511 

% 

37 

139 

29 1 

281 

15.9 

93 

100.0 

Sot>f*cis 

ft 

28 

70 

110 

85 

40 

20 

353 


79 

198 

31.2 

24 1 

11.3 

5.7 

1000 


The effect of multiple longitudinal obser¬ 
vations per individual subject within the dif¬ 
ferent age groups was evaluated to determine 
th^ possible dependence of multiple obser¬ 
vations in the same subjects over time. Only 
4.3®V of those with 3 longitudinal observa¬ 
tions were within the same age group. Those 
with 3 or more longitudinal observations had 
56*?t of such observations in 2 different age 
groups, 414V in 3* and 1 * 1 % in 4 different age 
groups. Thus, individual subjects tended to 
have observations spread throughout the age 
range. The number of observations was used 
as a covanate in later analyses to represent 
each subject’s contribution to those analyses. 

Method of Evaluating 
Functional Development 
The first approach chosen for the present 
evaluation of development (“growth") was to 
use all observations of all subjects over the 
age range in the longitudinal data set to find 
the best-fitting equations (using least squares). 
It is similar to the descriptive approaches used 
previously in this study (15-18). This approach 
is generally one of using unweighted obser¬ 
vations per individual to evaluate each pub 
monary function variable The best fitting 
composite equations were curvilinear relat¬ 
ing to age as well as to height. The robustness 
of the method was determined by comparing 
the values derived for subgroups with differ¬ 
ing numbers of observations. The values for 
those with only 2 points differed significantly 
from those with 3 or more observations. 
Therefore, those subjects with only 2 points 
were excluded from further analysis. The stan¬ 
dard error of the estimate for FEV, by num¬ 
ber of observations for 3 or more points 
ranged from 0.3 to0.5 L, and decreased with 
the number of observations; the multiple 
correlation coefficient increased with the 
number of observations from about 0.88 
to 0.95 for FEV,. Thus, predicted values fit 
the observed values very closely. The standard 
error of the estimate and the multiple corre¬ 
lation coefficient for the Vmax M did not 
have any trend related to the number of 
observations. 

Results of Composite Fitting 
Forms of the equations and resulting predic¬ 
tions derived from each of the other longitu¬ 
dinal subsets defined by the number of ob¬ 
servations were similar in appearance. Age 
(in decimalized years) and a height-height 
squared function (height in inches) were im¬ 


portant overall and in all age groups. In fe¬ 
males, height did not explain as much of vari¬ 
ance in parameters as it did in males. Age 
squared was important for the curvilinear 
model of Vmax^/FVC as well. 

The FEV, and FVC were greater for males 
than for females, but were a relatively con¬ 
stant function of height at each age. Thus, 
these variables and the Vmax**/FVC were 
evaluated using separate gender models (13, 
14) and/or using sex as a covariable (9). Com¬ 
pletely separate gender analyses did not con¬ 
tribute much useful information, and are not 
discussed at length. 

The final curvilinear equations that best 
explain the development of function between 
5 and 25 yr of age, based on 3 or more obser¬ 
vations per subject, were derived for FVC and 
FEV, (in liters), Vmax*, (liters/second), and 
Vmax w /FVC (in Uters/FVC seconds, as a 
measure of lung size-compensated flows). 
They are as follows: 

FVC - 9.17 + 0.054 Age - 0.373 Height 
+ 0.004 Height* (SEE * 0.495, R - 0.923, 

p < 0.001). 

FEV,,* 6.844 + 0.040Age - 0.281 Height 
+ 0.003 Height* (SEE - 0.431; R - 0.916, 

p < 0.001). 

Vmax*, - 5.489 + 0.056 Age - 0.221 
Height + 0.003 Heighri (SEE - 1.037, R * 
0.717; p < 0.001). 

Vmax»/FVC * 4.65 + 0.070 Age - 0.0019 
Age* - 0.116 Height + 0.0008 Height* (SEE 
- 0.296, R - 0.278, p < 0.001). 

Residuals 

Residuals from the equations did not corre¬ 
late with any further age, age-height, or size 
determinants (including sitting height and arm 
span). Average age did correlate with num¬ 
ber of observations; those hiving the largest 
number of observations were an average of 
I yr younger than those having only 3 obser¬ 
vations. This was not a factor that influenced 
the results. 

To analyze differences related to other fac¬ 
tors, residuals from the equations were de¬ 
rived for each observation and were expressed 
as proportions of the predicted values. Re¬ 
siduals are positive if above the predicted value 
and negative if below it; percent predicted can 
be calculated by adding 100*7* to the resid¬ 
ual. The residuals thus derived were normally 
distributed within age groups and overall, and 
age was not' a major factor in determining 
these distributions. 


Residuals for observations of ail subjects 
Had significant autocorrelations with at least 
3 prior observations of FEV, and Vmax* 
when entered hierarchically in multiple regres¬ 
sions. The simple correlation for FEV, values 
in adjacent surveys was 0.77, decreasing to 
Oj53 for values 3 surveys apart. Therefore, be¬ 
cause of autocorrelations, in most analyses 
we evaluated the effect of factors using only 
the end-point residuals for each individual 
subject; as representing that subject's outcome 
value. 

In addition, end-point residuals were caf 
culated for all subjects younger than 14 yr 
of age at entry with any spirometric test who 
were 12.4 to 26.4 yr of age at the time of their 
last spirometry (n * 440)j This allowed us 
to assess effects of early respiratory illnesses 
and effects of early exposures to tobacco 
smoke on end-point function. 

Statistical Comparisons 
Statistical analyses were performed on the 
DEC 10-Cyber 175 of the University Com¬ 
puter CemeT, using the SPSS statistical pack¬ 
age, the BMDP Package (19); and custom For f 
tran programs. 

Methods used included those of previous 
studies, such as multi factor analysis of vari¬ 
ance (ANOVA) and covariance (ANCOVA) 
(9), and more standard evaluations of sub¬ 
group differences (3, 5; 6, 8, 15-18). Age and 
sex were used as covariables, as in the analy¬ 
ses of Ware and coworkers (9), except when 
specified. Analyses evaluated the independent 
and interactive effects of various risk factors 
on lung function. As height is a function of 
age also, and as their effects on function are 
combined biologically between puberty and 
maximum lung growth (17), interactions of 
these 2 were used in these analyses. Because 
the method used is a variant of the repeated 
measures ANOVA (linear model) (19,20), the 
repeated measures ANOVA model was used 
as welli(19), in pan to verify the results from 
the method and the multifactonal ANOVA 
and ANCOVA. 

Individual Development Curve Model 
A second major approach assessing relation¬ 
ships between other factors and the multiple 
cumulative change in function used individual 
development (“growth^) functions of the term 
age times height* that best linearized the 
changes in lung function with time. This is 
the same term found to best describe the rate 
of change in adults in our study (18), and 
found previously with other descriptive meth¬ 
ods (17). The interaction is important bio¬ 
logically, as explained above. For subjects 
with at least 3 data points, standard errors 
of estimate derived in the process were used 
as measures of individual variability in anal¬ 
yses of the slope using the approach of Gold¬ 
stein (21, 22). In this approach^ developmen¬ 
tal (“growth rate") functions (or individual 
“slopes") are the dependent variables analyzed 
in generalized mixed linear models with 
covanablfcs (age, sex, height), risk factors of 
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interest, individual variability (as mentioned), 
numbers of observations (as mentioned), and 
years of follow-up. Hierarchical approaches 
were evaluated in these processes. One can 
then use variables representing change (in 
height, smoking, etc) as independent variables. 

Personal smoking and the standard airways 
obstructive disease (AOD) symptoms were 
derived from questionnaires completed over 
time bv the subject or by the parent when the 
subject was younger than 15 yT of age (12, 
14-18). 

Rasutta 

The three-dimensional solution of the 
equation for the FEV, is shown in figure 
1. It shows clearly the curvilinear rela¬ 
tionship of function with age and height 
(and curvilinearity with both together). 
The three-dimensional solution was of 
the same form for asymptomatic non- 
smokers (“normals”) and for others. At 
age 25, the results for these “other” sub¬ 
jects shows an apparent deficit in FEV, 
of 200 to 300 ml, approximately 6% of 
normal function. The equation is of the 
same form for the 2 genders, differing 
only in age contribution. The only vari¬ 
able showing a meaningful age-sex in¬ 
teraction was Vmax^/FVC. These data 
are shown in figure 2. 

Residuals of all pulmonary function 
variablfcs were significantly related to 


significant AOD symptoms in all age 14) was considered constant (always or 
groups; the most sensitive indicator was never), new, or remitted, as seen in table 
Vmax*> (tabic 2). Those who were asymp- 3. Ever diagnosed asthma and chronic 
tomatic had no significant variation of bronchitis had prevalence rates of 6.1 and 
residuals of lung function with age In respectively. About 45^» of those 

contrast, most symptom groups showed with either diagnosis had both. Those 
variability over the age range Those with who had physician-con finned asthma or 
cough appear to show an increasing loss chronic bronchitis throughout the study 
of function with age Asthma was related had significantly lower outcome FEV, 
to the worst lung function in each age and especially lower outcome Vmax*, 
group. than did those in the 3 other groups; the 

Outcome pulmonary function was trends in their function were not linear, 
evaluated by the longitudinal status of and these values after 20 yr of age were 
physician-confirmed disease, specifically still below nonnaL Those who developed 
asthma or chronic bronchitis. Longitu- asthma during the study had lower out- 
dinal status in those who entered the come Vmax*« (but not FEV\) than did 
study as children (initial ages zero to those who never had asthma or remit- 
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TABLE 2 

ALL RESIDUALS OF Vmax* BY SYMPTOMS BY AGE GROUP’ 




. 






AGE {yr) 



Toui Yp, 

(%| 

Evac H*i Sx 

AMt 

5.5-10 tl-15 

16-20 21-25 

Cough 

- 4 84 

♦ 0.51 -3,61 

-6 34 

-5 69 

239 

No 

♦ 1 43 





nafflui 

-7:20 

-1.93 -10.07 

-5 47 

-842 

158 

No 

+1:27 





WhM» 

-5.80 

-6.68 -5.11 

-5.81 

-6 12 

14 2 

No 

♦ 2 83 





Attack* tab 

-851 

-6Z2 -891 

- 10.44 

-6 17 

139 

No 

♦ 1 48 





Any Mtnmi 

-13.58 

- 17,97 -13.95 

-10.80 

15 69 

74 

No 

♦ 1.26 





* ANOVA YES war*u* NO. p < 0 001 tor Ym varauS No 

T Th* rvumb* of obaarvauon* (paraon* bma* Mm. tor paraona atofi 3 ♦ taatt) 

- 1.4*7 





TABLE 3 




OUTCOME FEV, AND ' 

< 

i 

| 

1 



CHILDREN < 14 YEARS OF AGE 





FEV, VmaAtt (unadjustad) 

Vmax* (adju*tad)T 


Numbar 

Maan SO Maan 

SO 


Maan 

MD aathma 






ANray* 

9 

-6.94 11.05 - 24 49 

19 20 


-24 8 

Naw 

24 

-0 69 16 74 -16.19 

2015 


-20.9 

ftamtsaioo 

17 

-0 77 16 04 -1.97 

29 71 


♦ 1.37 

Navar 

300 

♦ 3.02 13 37 *1.69 

25 13 


4 

OvaraJi (ANOVA) 


p » 0 07 

p < 0 0002 


p < 0.05 

Any asthma outcoma 


-2 36 p * 0 03 -17 43 

p < 0 0001 




* A ojuatao vatua* tor *91 mo Mtftma group* onty That* wort wgrvfc ar toy d#»arant 
and x qrxftcarsTy dr«*r*ot aga or owl o< w ewr *9 a 5 1 warm* 31 a 27 tor 
maana wrt >aiu»*d tor 09* by ANOVA. 09* or on** waa no* 1. • wgrutaant mart or 
Sax «m adfuatad tor. thougn not a agrHtcant awanabto. bacauaa or atq 


by aga or onaat 1 < S mut > 5) 
* and ‘ramaarona’ |p < 0 009) ao 
artact * ovaratt ANOVA vatu#* 


TABLE 4 

OUTCOME FEV, AND Vma^ RESIDUALS BY CHILDREN’S SMOKING AND 
SYMPTOMS CHILDREN < 14 YEARS Of AGE AT ENTRY (n - 237)* 






Smoking 






FEV, 



Vmax* 


Symp4oma 

Currant 

Ex 

Navar 

Currant 

Ex 

Navar 

Adjust adt 

Evar 

-2.77 

-1.59 

2 46 

-12.12 

-4.69 

0 10 

-433 


1.26 13.08 

p< 0.02 

2.71 

-1.90 8 93 

P < 0.04 

366 

2 43 

Adpjaiad* 

-149 

9.83 

296 

-7 44 

-0.04 

285 



I N rataaad 10 v "Wtotmry or w Aormaaon on a* vanattaa Ai p yatuaa by ANOVA 
t For FEV, adiuatad tor tympaon* paaon* amorung. paranaAOO and aga group 4 * • 

ta tgrufi ca w t) in ANOVA p < 0 077 tor adyuatad amotung; agxteao wracoona or amoving Mtn 
tom*(p - o.oaaj.ForVmax* adyuaiadmamaflaaatar aacn othar. aga aaa.paranta AOOano amo*- 
mg. and numbar or ooaarvabona (ai naj m ANOVA. p CO.031 tor a rt | wi fi amotang and p a 0 0*« tor 


ted. Those who remitted from asthma or 
chronic bronchitis during the study had 
end-point function between those who 
continued to have it and those who never 
had it. Outcome function in different 
asthma histories was si gni ficantly related 
to age of onset (in bivariate ANOVAs 
only), especially Vmax*, in those with “al¬ 
ways asthma”: mean outcome residuals 
were -31.2 for onsets <5 yr and -16.0 
for onsets > 5 yr. As average age of new 
onsets was significantly different also (6.9 
yr versus 3.1 yr for others), outcome 
Vmaxa* was adjusted for age (as seen in 
table 3). Sex and other factors were not 
significant covariates in effects of diag¬ 
nosis on function. 

After symptoms, smoking habits had 
an important effect on optcome. Resid¬ 
uals were significantly affected by smok¬ 
ing habits, both before and after adjust¬ 
ing for other factors (by muhifactor anal¬ 
yses of variance). The outcomes were 
related to symptoms within smoking 
groups, as shown for FEV, and Vmax*> 
(table 4); after adjustment for other fac¬ 
tors, symptoms were related still to out¬ 
come Vmax»«. Sex was a significant 
covariable in this multivariate relation¬ 
ship. The interaction between smoking 
and symptoms was significant for both 
outcomes, and those who toth smoked 
and had symptoms had the lowest end¬ 
-point residuals (table 4). The number of 
observations per individual subject, used 
as a covariable, was not significant. Pa¬ 
rental airway obstructive diseases (diag¬ 
nosed emphysema, chronic bronchitis, or 
asthma) had additional effects with chil¬ 
dren’s smoking and symptoms for end- 
point Vmaxjo, but not independently, the 
residual for current smokers with symp¬ 
toms and parental history was -14.41 
(95% confidence interval of -4.57 to 
-24.25). 

Parental smoking was not significant 
as a main effect within the multifactor 
analysis of variance for outcome FEV, 
or Vmax*>. However, it showed a signifi¬ 
cant interaction with subjects' smoking 
^children smoked if parents did) and pa¬ 
rental AOD in its effects on Vmax^fta- 
4>te 5). There were no such associations 
in subjects whose parents did not have 
AOD (not shown). Maternal smoking is 
significantly related to FVC (p < 0.001), 
producing a positive FVC residual: + 3.13 
versus -1.48 using sex-specific equations 
and +3.29 versus -1.36 using combined 
gender equations. This “inverse” relation- 
ship is found in all age groups, but is sig¬ 
nificant primarily in the older age groups 
and in females. 


For size-compensated flows (Vmax**/ 
FVC), females had higher predicted val¬ 
ues between ages 13 and 25 (figure 2). 
Size compensation removed any further 
effects of age on residuals. Subjects’ 
symptoms were still highly correlated 
with the outcome residuals (table 6), 
regardless of gender. There was a signifi¬ 
cant relationship for subjects’ smoking 
with outcome residuals, using sex-specific 


or combined models. (The effect of ever 
smoking was seen more in those 15 to 
25 yrof age at end point.) Parental smok¬ 
ing was not significant, though^moking 
habits of mothers showed a tfena; (Ma¬ 
ternal smoking was significant indepen¬ 
dently (p < 0.012) only when all observa¬ 
tions were evaluated together!) There were 
no significant interactions in the mul¬ 
tifactor ANOVA. 


t 
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TABLE 5 


OUTCOME Vmu* RESIDUALS BY CHILDREN'S SMOKING*. 
PARENTS SMOKING. AND PARENTS' AOD* CHILDREN < 14 YEARS 
OF AGE AT ENTRY* 



Total* 


Patents with AOD 

Smoking 

Currant 

Now 

Current 

IwVff 

Pargnts Smoking 

Mofhar* 

- 9 131 

44* 

-9 09 

236 

Fathar only 

-5 73 

1.28 

-3.23 

-11.04 

Fathar* 

-7 07 

071 

-6 53 

-700 

Naitnar 

-0 40 

5.86 

0Z5 

1368 

ANOVA (vs nanftar) 

Method and tatnar or#y ns. 
F«twrt, p 4 0.03 

ns 

< 0 06 



* b-tmour* not 

^ As AlQNIMO slWpftySSfflS. CAWW bfWCMS. Of SStlVna 

* Agn 12 * » »*yr « and port 

1 Muc*d N rMMd to tvtktoitoy ot o WtiMWO on ON vtnMto* and 090 (170) 
*£/l tgoust mwkmg 

I Mto conkosocs **»***< « -1*5 to • 1 ** 1 . — oAw onoowpoos 0 


TABLE 6 


OUTCOME Vm*x*/FVC RESIDUALS RELATED TO CHILDREN'S SMOKING AND 
SYMPTOMS. AND PARENTAL SMOKING (n - 389): 

CHILDREN < 14 YEARS OP AGE AT ENTRY 


Adjusted todapandant Factors* 

Adjust *d 

Ma*n* 

p VMuas 

Sudfocts symptoms 

Y«s 

-3 24 


NO 

-0.67 

< 0.014 

Subjects smoking 

Evr 

- 10.02 


Nsvsr 

- 1i25 

< 0.013 

Mqdters smokmgt 

Currant 

-4.41 


Ex 

-3.61 


Navar 

-0 40 

> 0 42 

OvaraH ANOVA 


< 0 . 001 


* Yitost m*»n sfteets sdfustod tor an (p < O.OOi). ago (not. and sscn ochor tomsfs omoamff rsm ovso 
(aa Mcona«iam tmnd) 
t A rl t u ao d tor tsthsr* smoking aa won 


TABLE 7 


INDIVIDUAL GROWTH RATES* OF FEV, AND Vmu* BY RISK FACTORS* 
COVARIABLES * AND INDIVIDUAL VARIABILITY* 



NumterV 

Adjusted FEV, 

Adjusted Virtu* 

Smoking 

Currant 

52 

0.74 

0.65 

Ex 

20 

0.95 

092 

nm. 

137 

0*3 

p • 0.055 

0.93 

p - 0.005 

Symptoms 

Yaa 

(79) 

- 

0.74 

NO 

(130) 

ns 

0.93 

p • 0.014 


* Coofecto* M (ags a HI*) * 1*00. 

* Bmotong and s ymptoms adpuatod tor aacn mm aa «*■ aa to* mm factors (toctodmg psronta* 

•nwtong. wtoeh was not stgnrk ta mj 

* Aqs and •*■ ( aiy a Ac a np. 

* Utong SEE f r crroiswd wtto # too* and aga. not wito raft tocaon). and nuwtoor 0 1 nmin town*. 

(Conototod w«to aga. yMn ot to*ow up. and w*» wadjiiart FEV,.) 

1 RsOucrt o o c o u oa ot toco ol data on ndmduat amotong to g., ogaa undo* II yr) or s ymptoms 

The use of the repeated measures fitted developmental curves analyzed by 
ANOVA Unear model confirmed the multifactorial ANOVA/ANCOVA -fir 
findings posited. Thus, residuals from the the more general Unear model. 


Results of Individual 
Development Mode! 

As previously stated, individual “growth” 
rates (“slopes”) were evaluated following 
the approach of Goldstein (21, 22), us¬ 
ing multifactorial analysis of covariance 
Standard errors of estimate (SEE) were 
used as a covariate for individual vari¬ 
ability (within and between surveys). The 
SEE were significantly correlated with 
the number of tests and age but not with 
any of the risk factors. The number of 
observations per individual was used as 
a covariable also. It was correlated with 
age (and, of course years of follow-up). 
It was correlated significantly with un¬ 
adjusted mean FEV, in the different sur¬ 
veys, though without trend, indicating 
random survey differences. The in¬ 
dividual “slopes” were not correlated with 
initial values, but were correlated with 
end-point (final follow-up) values. Af¬ 
ter each explanatory variable was adjust 
ed for covariables and ^^.explanatory 
variables, only symptom!?and smoking 
were highly related to individual slopes, 
especially those of Vmax M (table 7). 

Discussion 

The composite unweighted method was 
used to obtain the best descriptive fit for 
lung function development (ie., growth) 
curves. It was not developed only to look 
for deviations/variations from time 
trends, as with time series models. The 
method used was found to be very ro¬ 
bust. It adjusted for age and for body 
size (standing height k after which other 
measurements (sitting height, arm span) 
did not contribute to explaining func¬ 
tional measurements. Unfortunately, we 
did not make chest measurements, which 
have been found to help explain the 
growth curve in late adolescence and 
early twenties (8). The method is not in¬ 
tended as a reference formula. It is the 
mathematical best fit to the data and is 
not assumed to adequately describe bio¬ 
logic growth, although the solutions par¬ 
allel our previous descriptive results (18V 
The method provided very similar 
results for those with different numbers 
of observations, as long as subjects have 
3 or more observations over an 8- to 12- 
yr period. Individuals with only 2 points 
do not provide sufficient information 
and do not fit the curvilinear solutions. 
Over lengthy time periods, observations 
on an individual subject appear to be dis¬ 
tributed over the entire age range, thus 
minimizing dependency of observations. 
(Thus, there were no major losses of 
degrees of freedom or increased inter- 
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dependence of the observations based on 
the number of values used for an in¬ 
dividual subject.) We used values that 
could be examined at the follow-up end 
point and at onset. The use of residuals 
was very convenient statistically in that 
it led to variables that were approximately 
normally distributed, around a zero av^ 
erage, and was equivalent to looking at 
percentages above normal and percen¬ 
tages below normal in evaluating other 
factors affecting the development curve. 

The use of the repeated measures 
ANOVA linear model (19, 20) as an ad¬ 
junct to the composite-residual method 
showed that the descriptive fit, used for 
physiologic purposes/results yielded a 
good fit to a general linear model. The 
use of ANOVA/ANCOVA to analyze 
differences, as per Ware and coworkers 
(9), was a very robust and pragmatic ap¬ 
proach. Results that can be provided as 
actual decrements of function are far 
more understandable than are odds ra¬ 
tios or chi-square results. 

. Residuals were used to examine the 
phenomena of “tracking” the subjects* 
values over time. Thus, the relation of 
each one to “average” growth in this 
population was evaluated by looking at 
any change in the relation of the devia¬ 
tions from average (i.e., these residuals) 
over time. Normal subjects showed neg¬ 
ligible mean change in their residuals 
from initial to end point per annum. The 
grand mean change in residuals was only 
0.26%. even including those with symp¬ 
toms, smoking, and other risk factors. 
Only active smoking disturbed the track¬ 
ing with any significance (p * 0.099), 
with current smokers having an adjusted 
mean change of -0.64 ^Symptoms, 
often present initially, only produced 
a -0.19% change (p * 0.17). Other 
factors had no apparent influence on 
tracking. 

The number of observations related 
sign: cantlybut without pattern only to 
survey. In all methods of analysis, the 
number of observations was used as a 
covariate as a measure of individual varir 
ability; its lack of significant contribu¬ 
tions was considered another indication 
of the strength of this approach to model¬ 
ing. The SEE derived for each subject 
with > 3 values represents a measure 
of within-individual variability. It was 
found to be unrelated to any risk vari¬ 
ables, and related only to age and the 
number of observations. The SEE was 
used as a covariate in analyses of growth 
rates (individual “slopes”), using a lin¬ 
ear model after the fashion of Goldstein 


(21, 22), to determine contributions of because parents do. As found by many 
risk variables to changes over time. The others, there was a significant! relation 
results therefrom substantiated previous ship between the smoking habits of,pa r . 
results, as wdlias providing another meas- ents and those of their children, especially 
ure of outcome. between those of the same gender. Like- 

Size-compensated flows (expressing wise, mpthers and fathers who smoked 
Vmax w in FVC seconds) showed the same usually had spouses who smoked, 
pattern with age and by sex (figure 2) as further, parental smoking had an in. 
shown in infants and smaller children (10, fluence^on the ^children’s pulmonary 
23); females had higher flows for given functiomoutcome. It was most evident 
volumes than did males after puberty, in current smokers, symptomatic sub¬ 
even though males had higher volumes, jects, and those with a parental history 
For size-compensated flows, sex-specific of AOD. Mothers’ smoking was signifi- 
analysis yielded the same results as those cant, as found by others (9, 28)* but we 
obtained for the combined group using also found fathers’smoking to beimpor- 
gender as a covariable Flows per se tantL Parental smoking did not influence 
(Vmaxso) did not differ significantly by height at any age as suggested by these 
gender As discussed by us previously (10, previous studies. We found that subjects 
14-18, 23) and by others (3, 8, UK the did have increased volumes (FVC) and 
FEV, and flows do measure airway and decreased flow^if parents smoked. The 
parenchymal changes, and one can dis- effect did not differ by gender and was 
criminate differential contributions, es- noted even in 5 to 7-yr-olds, although the 
pecially when evaluating flow as a ratio magnitude was greater in older children 
to vital capacity. and young adults. It is interesting that 

Symptomatology contributed as well preliminary studies on infants (29) sug- 
to outcome function as another predic- gest that children of smoking mothers, 
tor of risk, but not to individual rates especially males, may have elevated func- 
of growth. Outcome function, of course, tional residual capacities even shortly af- 
was very well correlated with outcome ter birth. We have observed also that 
symptomatology as well, and to individ- adult male smokers who show a rapid 
ual growth rates. Outcomes related spe- subsequent functional decline are likely 
cifically to incidence and remission of to have well-preserved FVC values at an 
asthma and chronic bronchitis. As shown early stage of their illness (30). The mech- 
by others (11), abnormalities of flow anism underlyingthes^obsei^ttiotis re- 
(i^., the Vmax w ) persist longer, into mains unclear, but it would appear that 
adulthood. increase in lung volume could be an 

Further children’s symptomatology warty manifestation of the effects of ac- 
(by parental response and self response) live or even passive*smoking, related to 
was related to parental history of AOD. stimulation effects of nicotine, or growth 
The latter did relate to outcome of symp- compensatory effects related to CO ex- 
tomatology, and it modified effects of posure and/or possible loss oMung 
other factors (see below). Initial meas- elasticity. 

ured pulmonary function was influenced Our robust method of lung function 
by both parental factors and the chil- growth showed effects on function of var- 
dren’s respiratory history, and growth in ious risk factors that are highly compat- , 
function was affected by these factors, ible with those seen using the traditional 
the level of their initial function, and later general mixed linear model and addi- 
symptomatology. This information con- tional linear model of individual growth 
firms longitudinally, in a preliminary (21, 22). Thus, the linear models used are 
fashion at least, that our hypotheses that consistent in this population in reflect- 
childhood respiratory trouble, effects of ing actual growth of pulmonary function 
key respiratory illness, and the effects of and the factors that affect it. 
parental/familial factors are all impor- In conclusion, pulmonary function 
untin growth of pulmonary function in growth is significantly related to chil- 
children (24-26). dren’s respiratory symptoms and disease. 

Smoking had a major effect on growth and their smoking; parental/ familial fac- 
of pulmonary function as well, aadjjtere tors are important in some subgroups as 
was an interaction of smoking habttfof well Surprisingly, symptoms/disease are 
{$>e subjects and symptomatology on out- not integral in analyses in other studies 
come of pulmonary function. Ex-smok- (9, 28), nor are the interactions of int¬ 
ers may have shown a rebound phenome- portant risk factors. The effects of pa- 
no n, seen in analyses of young adults rental smoking are somewhat inconsis- 
(27). Children appear to smoke in pan tent; in some cases a maternal smoking 
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cffcci is noted only within other risk fac¬ 
tor subgroups. The independent and in¬ 
teractive effects are seen as perturbations 
in the growth curve as well as in out¬ 
come function (FVC, FEV*, Vmax w . 
Vmax»o/FVC). Flows (Vmax M ) and size- 
compensated flows (VmaWFVC) were 
the most sensitive, and often showed per¬ 
sistence of decrement related to early 
events, 
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